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The 7-ray emission of the blazar Markarian 421 above 250 GeV has been observed by the CAT Cherenkov 
imaging telescope since December, 1996. We report here results on the source variability up to April, 1998, 
with emphasis on the 1998 campaign. For the flaring periods of this year, the energy spectrum was derived 
from 330 GeV up to 5.2 TeV: it is very well represented by a simple power law, with a differential spectral 
index of 2. 96 ±0.13. 
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Total time ON source = 73.9 hours 

Signal = 2877 ± 380 
7 / min = 0.65 ± 0.09 
a = 7.6 
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1 Introduction 

Markarian 421 (Mrk 421) is one of the three extragalactic sources detected to date from the North- 
ern hemisphere at TeV energies, along 
with Markarian 501 and 1ES2344. In 
particular, it is the closest BL Lac ob- 
ject (z = 0.030), as well as the first 
one discovered in the very high energy 
(VHE) range: shortly after its detec- 
tion between 50 MeV and a few GeV 
by the EGRET instrument on board 
the Compton Gamma Ray Observa- 
tory (Lin, et al. 1992), the Whipple 
group observed it at a 6. 3a level at 
TeV energies (Punch, et al. 1992). 
Since then, it has remained one of the 
most studied blazars, and was the ob- 
ject of several multi-wavelength ob- 
servation campaigns (see, e.g., Ma- Pointing Angle, a ( ) 
comb, et al. 1995, and Takahashi, et Figure 1: Distribution of the pointing angle a (see Sect.]2\for definition) of Mrk 421 
al 1996) for all data between December 1996 and April 1998. The signal is clearly seen in the 

TlfSt r)lYL 

Since it started operation in Autumn 

1996, the CAT (Cherenkov Array at Themis) telescope has been observing Mrk 421 regularly. The next sec- 
tion describes briefly the detector and the analysis method used to extract the signal. Then, the Mrk 421 
data sample and the corresponding light curve up to April 1998 are presented in Sect. |3] Finally, the spectral 
properties observed in 1998 are examined in Sect.|4] and a discussion is given in Sect. [5] 

2 The CAT detector: characteristics and analysis method 

The most complete descriptions of the CAT experiment and of its specific analysis method can be found 
in (Barrau, et al. 1998) and (Le Bohec, et al. 1998), respectively. We just recall here the few characteristics 
which are necessary for the understanding of the following sections. 

Located on the site of the former solar plant Themis (French Pyrenees), the CAT detector uses the now- 



validated technique of Cherenkov imag- 
ing. The Cherenkov light emitted by 
the secondary particles produced dur- 
ing the development of atmospheric 
showers is collected through a large 
mirror which forms their image in 
its focal plane. The CAT camera, 
placed 6 m from its 4.7 m diame- 
ter mirror, has a 4.8° full field of 
view, which is comprised of a cen- 
tral region of 546 0.12° angular di- 
ameter phototubes in a hexagonal ma- 
trix and of 54 surrounding tubes in 
two "guard rings". Its very high def- 
inition, combined with fast electron- 
ics, is an efficient solution of the two 
main problems one has to face with 
such 7-ray detectors: firstly, the detec- 
tion threshold, which is determined by 
the night-sky background, is relatively 
low (250 GeV at Zenith); secondly, 
the capability of rejecting the enor- 
mous cosmic -ray background (protons, 
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Figure 2: Mrk 421 nightly integral flux levels above 250 GeV, for observations 
nuclei) is improved by an accurate between December 1996 and April 1998. The spectral shape as derived in Sect&ias 
analysis based on the comparison of been assumed to estimate the integral flux for observations far from the Zenith. The 

individual images with theoretical vs\&sn^ as ^ ie ^ ^ ne r epresents the mean flux over the two years. 

7-ray images, i.e. by the fit of an analytical model which uses the whole information contained in the images' 
longitudinal light profile. This fit allows the direction and energy of 7-ray events to be determined with good 
accuracy, with an angular resolution per event of the order of the pixel size, and an energy resolution about 
20% (independent of energy) when selecting showers with an impact parameter less than 130 m. 
Two main cuts are needed for the 7-ray selection. After keeping the highest-energy triggers (total charge 
Qtot > 30 photo-electrons), a first cut on the probability given by the fit P{x 2 ) > 0.35 is applied, which 
selects the images with a "7-like" shape. Since 7-ray images are expected to point towards the source position 
in the camera, a second cut a < 6° is used, where the pointing angle a is defined as the angle at the image 
barycentre between the actual source position and the reconstructed angular origin of the image. As a result, 
this procedure rejects 99.5% of hadronic events while keeping 40% of 7-ray events. A source like the Crab 
nebula, which is generally considered as the standard candle for 7-ray astronomy, can be detected at a 4.5<r 
level in one hour and localised within 1 to 2 arcmin. 

3 Data sample 

The 7-ray emission of Mrk 421 above 250 GeV has been observed by CAT from December 1996 onwards. 
Here we consider data taken up to April 1998. A selection based on criteria requiring clear weather and stable 
detector operation has been applied. This leaves a total of 73.9 hours of on-source (ON) data, together with 
33 hours on control (OFF) regions, for a zenith angle band extending up to 40°. Fig.^shows the distribution 
of the pointing angle a for the whole sample. With the cuts quoted in Sect. 13 the mean 7-rate is found to be 
of the order of one third of that of the Crab nebula, with a total significance reaching the 7.6a level. After 
Mrk 501, this makes Mrk 421 the second extragalactic source detected by CAT. 

As can be seen on Fig. |2j the emission of Mrk 421 changed significantly between 1996-97 and 



1997-98: almost quiet in the first period (with 
a mean integral flux (/>(> 250 GeV) = 1.94 ± 
0.65 x 10~ n cm~ 2 s -1 ), the source showed small 
bursts in the second period, together with a higher 
mean activity (</>(> 250 GeV) = 6.05 ± 0.54 x 
10~ n cm~ 2 s _1 ). Fig. |3]is a zoom on the small 
flare of the night between 24 and 25 March, 1998, 
which reached a maximum of ~ 2 times the 
Crab emission level. Three days after this flare, 
a VLBI observation showed a strong decrease of 
the radio emission polarization and intensity from 
the core of the galactic nucleus (Chariot, et al. 
1999), in comparison with a first measurement 
performed three weeks earlier. The insufficient 
time-coverage between the radio and TeV obser- 
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Figure 3: Light curve of the night between 24 and 25 March, 1998 
(MJD 50896/50897). Each point represents a run of ~ 30min. 



The data used in this section 
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vations unfortunately precludes any further interpretation. 

4 Mrk 421 spectrum in 1998 

A spectrum was derived for the flaring periods of Mrk 421 in 1998 
have been further limited to zenith an- 
gles < 16°, i.e. to a configuration 
for which the detector calibration has 
been fully completed. By minimizing 
systematic effects, this allows a ro- 
bust spectral determination. Another 
favourable factor is the low night-sky 
background in the field of view, due 
to the lack of bright stars around the 
source. 

The spectrum was derived in two steps: 
firstly, the number of events for each 
energy and zenith angle bin within 
the cuts quoted in Sect. |2] was de- 
termined for ON and OFF runs; sec- ' Enemy (TeV> 

ondly, taking into account the trig- 

, , rr- • ■ ,, ., Figure 4: Differential flux of Mrk 421 for the flaring periods in 1998. Individual 

ger and cut efficiencies as well as the ... , . , . ,. . , , J , . . , „ r „ ,. 

p intensities per bin are only indicative and the error bars are statistical. 1 he full line 

energy resolution OI tne telescope, a re p resen [ S the power law parameters estimated by a likelihood method, in which only 

maximum likelihood estimation of the the fullfilled circles were used (see text). This method gives also the 2a error "box", 

Spectral parameters has been applied whose shape reflects the correlation between the parameters. This latter vanishes at 







for the simple power law hypothesis 



the decorrelation energy E<j = 580 GeV. 



(Ho) d^/d-E/TeV = ^o-^ToV- ^ S can ^ e seen i n Fig- IH tn i s accounts very well for the observed spectrum 
between 330 GeV and 5.2 TeV, which was the energy range considered. With a correlation coefficient of 
—0.68, the fitted parameters are: 

0o = 1-96 ± 0.20 stat ± 0.12 s y s " MC ±o.59 SyS ~ atm xl0~ n cm^s^TeV- 1 and 

7 = 2.96 ± 0.13 stat ± 0.05 sys " MC . 

Systematics errors due to limited Monte Carlo statistics (noted "sys-MC") in the determination of the effectice 

detection area affect both cpo and 7. The second systematics errors (noted "sys-atm") come from the uncer- 



tainty on the absolute energy scale. They are due to variations of the transparency of the atmosphere and affect 
only the absolute flux. 

5 Discussion 

The spectral shape of Mrk 421 observed by CAT in 1998 is fully compatible with recent results from 
the HEGRA collaboration (Aharonian, et al. 1999). It is however in contrast with the former behaviour of 
the source at the time of the 1995 and 1996 flaring periods, as observed by the Whipple group, who found 
7 = 2.54 ± 0.03 stat ± 0.10 sys (Krennrich, et al. 1999). The higher value of the differential spectral index 
found here thus seems to favour a correlation between the intensity level and the spectral shape. Actually, such 
a correlation has been recently reported by CAT on the extreme blazar Mrk 501 (Djannati-Ata'i, et al. 1999; 
see also Tavernet, et al. 1999). Since the spectrum of the latter is now known to be curved for very intense 

emissions, the hypothesis {Ti\) of a curved shape d^/dE'xeV = ^o^TcV ° gl ° T was a ^ so cons idered 
for Mrk 421. The fitted curvature term is j3 = 0.28 ± 0.49, compatible with zero. The likelihood ratio A = 
—2 x log c(n °l , which gives an estimate of the relevance of H\ with respect to Ho and behaves asymptotically 

like a \ 2 with one d.o.f., is 0.34, corresponding to a chance probability of 0.56. This confirms the absence of 
any obvious spectral curvature for Mrk 421, as already indicated by Krennrich, et al. 1999. In the framework 
of leptonic models (Ghisellini, et al 1998), which succesfully explain the Mrk 501 spectral energy distribution 
in the X-ray and VHE 7-ray ranges (Djannati-Ata'i, et al. 1999), the present result implies that the peak energy 
of the inverse Compton contribution of Mrk 421 is significantly lower than the CAT detection threshold. This 
is not surprising since the corresponding synchrotron peak is lower than that of Mrk 501, and since leptonic 
models predict a strong correlation between X-rays and 7-rays. Moreover, such a correlation was directly 
observed on Mrk 421 in Spring 1998, during a coordinated observation campaign involving ground-based 
Cherenkov imaging telescopes (Whipple, HEGRA, and CAT) and the ASCA X-ray satellite (Takahashi, et al. 
1999). 
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